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Abstract 

Background: The intestinal microbiota is a complex and diverse ecosystem that plays a significant role in 
maintaining the health and well-being of the mammalian host. During the last decade focus has increased on the 
importance of intestinal bacteria. Several molecular methods can be applied to describe the composition of the 
microbiota. This study used a new approach, the Gut Microbiotassay: an assembly of 24 primer sets targeting the 
main phyla and taxonomically related subgroups of the intestinal microbiota, to be used with the high-throughput 
qPCR chip 'Access Array 48.48', AA48.48, (Fluidigm®) followed by next generation sequencing. Primers were designed if 
necessary and all primer sets were screened against DNA extracted from pure cultures of 15 representative bacterial 
species. Subsequently the setup was tested on DNA extracted from small and large intestinal content from piglets with 
and without diarrhoea. The PCR amplicons from the 2304 reaction chambers were harvested from the AA48.48, 
purified, and sequenced using 454-technology. 

Results: The Gut Microbiotassay was able to detect significant differences in the quantity and composition of 
the microbiota according to gut sections and diarrhoeic status. 454-sequencing confirmed the specificity of the 
primer sets. Diarrhoea was associated with a reduced number of members from the genus Streptococcus, and in 
particular S. oloctolyticus. 

Conclusion: The Gut Microbiotassay provides fast and affordable high-throughput quantification of the bacterial 
composition in many samples and enables further descriptive taxonomic information if combined with 454-sequencing. 
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Background 

Immediately after birth the mammalian gastrointestinal 
tract is colonized by a complex and diverse microbial eco- 
system. The bacterial invasion and the following gut micro- 
bial composition has an enormous impact on its hosts 
health and well-being [1]. To gain a better understanding of 
this complex ecosystem, culture-independent methods are 
essential, as a considerable fraction of the intestinal micro- 
biota has not yet been cultured [2]. One of the ongoing 
controversies is how to study the bacterial composition 
in complex ecosystems. To date, some of the widely used 
approaches to characterize the intestinal microbiota are: 
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metagenomics, phylogenetic microarrays, DNA fingerprint- 
ing techniques, and qPCR [2-5]. These methods provide 
different degrees of taxonomic as well as quantitative infor- 
mation on the microbiota [5]. Nonetheless, variation in 
technical procedures and differences in data treatment and 
interpretation makes it challenging to compare results be- 
tween studies. Also, the expenses and time consumption in 
relation to labour intensity and data analysis vary greatly. 
Especially metagenomic approaches are receiving increased 
attention in the study of microbial communities as a result 
of their shorter sequencing speed, extended read length, 
and lower costs [5,6]. However, the enormous amount of 
data generated becomes cumbersome to analyse, and re- 
quires lots of dedicated time as well as expertise to manage 
[6]. The Access Array 48.48, AA48.48, (Fluidigm Corpor- 
ation, South San Francisco, CA, USA) creates an affordable 
link between high-throughput qPCR and next generation 
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sequencing (NGS) and provides manageable data with valu- 
able quantitative and taxonomic information. 

Since the 1990s, qPCR has been applied widely due to 
its quantitative precision, high specificity and sensitivity, 
broad dynamic range, good reproducibility, and rela- 
tively low costs [7]. In general, qPCR is quick to perform 
with a low to medium throughput, since most qPCR 
platforms have a capability format of 96 or 384 [7]. The 
high-throughput qPCR chip AA48.48 combines 48 detector 
inlets and 48 sample inlets by interconnecting channels into 
2304 individual reaction chambers (singleplex) http://www 
fluidigm.com/access-array-system.html In contrast to most 
qPCR platforms which require reaction volumes between 5 
to 100 ul [7], the AA48.48 operates with a reaction volume 
of approximately 35 nl [8]. By tagging primers at their 
5 '-end (using the 'The Access Array 4-primer amplicon 
tagging strategy [8]), and including a 454 Barcode Library 
(454BL) with the samples, amplicon generation and library 
preparation is achieved in the same reaction. Afterwards, 
amplicons can be harvested directly from the AA48.48 
sample inlets where the respective samples were initially 
loaded. The unique barcodes and specific primers enable 
segregation of the pooled samples in NGS analysis later 
on. The AA48.48 process time from start to finish is 
approximately five hours. After each qPCR run Fluidigm 
Real-Time PCR Analysis software (Fluidigm Corporation) 
generates a heatmap constructed as a grid of 48 x 48 small 
squares presenting all the 2304 primer sample combina- 
tions. The heatmap depicts the Cq-value for each reaction 
according to a colour scale, providing a good overview of 
each sample and the possibility to compare the bacterial 
profiles across all the samples. 

This study describes the design, optimization, and valid- 
ation of the Gut Microbiotassay: a primer setup consisting 
of 24 primer sets targeting the main bacterial phyla of the 
intestinal microbiota at various taxonomic levels, to be used 
with the AA48.48 in combination with NGS. Furthermore, 
it demonstrates the applicability of the Gut Microbiotassay 
on luminal content collected from the small and large 
intestine of piglets of different diarrhoeic status. Finally, 
it validates the specificity of the Gut Microbiotassay 
by sequencing these amplicons. To our knowledge, the 
AA48.48 has not previously been used to investigate 
microbial communities. 

Methods 

Development of the Gut Microbiotassay 
Primer design 

Inspired by Rajilic-Stojanovic et al. [9], the primer setup was 
designed to target the ribosomal RNA genes (16S or 23S) 
of the major bacterial groups present in the mammalian 
intestinal microbiota, including the phyla: Firmicutes, 
Bacteroidetes, Proteobacteria, and Actinobacteria [10]. 
To gain insight into the bacterial composition, different 



taxonomic levels were represented: domain, phylum, class, 
family, genus, and species-specific primers (Table 1). 

Primer specificity was checked in silico with the probe 
match' facility of the Ribosomal Database Project, release 10 
(http://rdp.cme.msu.edu/) [11] and the check' application 
of ProbeCheck (http://131.130.66.200/cgi-bin/probecheck/ 
content.pl?id=home) [12] using the programs' default set- 
tings. Judged by the search results from these online da- 
tabases, primers were renamed according to their main 
target. Some of the published primers were slightly modi- 
fied to improve their specificity. New primers were designed 
using the 'Probe Design Tool' of the ARB Software package 
(http://www.arb-home.de/) [13], or Primrose [14] using de- 
fault options, and further validated based on BLAST search 
(NCBI). Primers were designed to have approximately same 
length of nucleotides, GC-content, a minimum number 
of degenerate bases, and to produce amplicons between 
200-500 bp long, compatible with the read length of the 
454 GS FLX Titanium Sequencer. Similar properties were 
important, given that all the primers had to function under 
the same conditions when run on an AA48.48. 

All primers, except for the species-specific ones, 
were tagged at their 5 '-end to enable incorporation of 
a 454BL necessary for 454-sequencing: forward tag: 5'- 
ACACTGACGACATGGTTCTACA-3 ', and reverse tag: 
5 -TACGGTAGCAGAGACTTGGTCT-3 ' in accordance 
to the Access Array™ System User Guide' [8]. Primers 
were purchased from Eurofins MWG Synthesis GmbH 
(Ebersberg, Germany) and stored at -20°C. 

Empirical testing of primers 

All primers were tested against 15 strains of pure-cultured 
reference bacteria, representing targets for one or more 
of the different primer sets (Table 1). Thus the reference 
bacterial DNA functioned as both positive and negative 
controls for the individual primer sets. 

DNA extraction Chromosomal DNA from cultured refer- 
ence bacteria (Additional file 1: Table SI) was isolated using 
the Easy-DNA™ Kit (Invitrogen, Carlsbad, CA, USA), 
for further details on the DNA extraction protocol see 
Additional file 2. DNA was precipitated with ethanol, and 
resuspended in 60 ul TE buffer. The reference bacteria were 
cultured as recommended by DSMZ (www.dsmz.de/). 

An interplate calibrator (IPC) was included in all 
AA48.48, consisting of bacterial DNA extracted from ~ 
100 mg colonic content from a healthy conventional pig, 
14-week-old, Danish landrace. Intestinal content was col- 
lected immediately after euthanization and frozen at - 80°C. 
A 10% PBS suspension was made from the intestinal 
content, and from here on, the protocol was identical 
to the one used for the reference bacteria. 

DNA concentration and purity were assessed by 
the 260/280 nm-ratio using the Nanodrop" ND-1000 
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Table 1 Target-specific primer sequences constituting the Gut Micro biotassay in order to target main bacterial phyla 
and bacterial groups and species of high interest in the mammalian intestine 

Main target rRNA Sequence (5'— > 3') E. coli Size, Reference Reference 







gene 




position 


bp 




bacterium 


Domain Bacteria A V2-V3 


Forward 


16S 


AGAGmGATCCTGGCTCAG 


7 


336 


Liu et al. [30] 


1,2,^,^,6,7,8,9, 




Reverse 


16S 


CTGCTGCCTYCCGTA 


342 




Liu et al. [30] 


10,11,13,14,15 


nnmsin R?irtpri3 R \/A-\/^ 
U'Ullldlll DdLlclld D vH vD 


rui vvui u 


1 6S 


CAGCAGCCGCGGTAATAC 

LnuLnULLOLuu I AA I Av_ 


518 


389 


OLI IVVIcytrl (y. IcUUc ID I J 


1 ,Z,D,H,D,V,/ ,0,y, 




Reverse 


16S 


CCGTCAATOCmGAGm 


906 




Schwieger & Tebbe [31] 


10,11,13,14,15 


Phvliim Firmin itpc. 
riiyiuiii i ii 1 1 iilulcj 


Fnr\A/n rrl 

1 Ul VVU 1 u 


1 6S 


CTGATGGAGCAACGCCGCGT 

^_ I \Jr\ I VD^AVD^.. AA^.VJ^.^.*J^.*J I 


385 


429 


nddixci ibci i ct ui. \_~>/-\ 


67011 1114 




Reverse 


1 6S 


A C A C YT A G Y A CJC ATCC-iTTT 


81 3 




K/li"ihlinn pf nl R31 

IVIUIIIIIiy Ct Ul. 




flacc Rprilli 
^Jdbb DdLIMI 


Fni r\ a in rn 1 
rui VVUI U 


1 6^ 

I Dj 


c-,c a ct accc a atcttcccc 


jjj 


461 


FpIcI/p pf nl 
rfcrlbKc ct Ui. ID H] 


71114 




Reverse 


16S 


ACACTOGCACTCATCGm 


813 




Modified from 
Muhling et al. [33] 




ucl lUb 01 ILcI ULULLUb 


Fni r\A/n rn 1 
rui VVUI U 


1 6^ 


CCCTA ACCTRCCC ATCAC A 


1 9^ 
I Zj 


J/J 


K/lnHifipH frnm Rphr Pt nl [3^1 
IVIUUIIIcU IIUIII Del II ct Ul. IDDj 


J 




Reverse* 


16S 16S 


GTOCTCTCATCCTOTO 
ACCGTCAGGGGACGTOAG 


449 466 


342 


Modified from Behr et al. [35] 
Modified from Behr et al. [35] 




Genus Lactobacillus 


Forward 
Reverse 


23S 
23S 


GCGGTGAAATOCAAACG 
CCC A CCTT A A CTCCTC AT 


774 
989 


216 


This study, Linux ARB [13] 

Thk cti irlv I iniiY ARR H 31 

1 1 lib blUUy, Lll IUA AnD [ 1 OJ 


11 


otrl lUb DLI cfJLULULLUb 


Fni r\ a in rn 1 
rui VVUI U 


1 ^ 


fTW A CC A C A A A CCC A CC CCT 


too 


DD / 


Thk cti iHv riiKfpil\A/9 T361 

II lib bLUUy, ^.lUbLdlVVZ IDOj 


14 




Reverse 


1 6S 


AACCRYCYAACACCTACC 

nnuun i ^_ i aa^.a^.^. i a^o^.. 


824 




Thk cti irlv riiKt?il\A/9 T3A1 

1 1 lib bLUUy, ^.lUbLdlvVZ iDUj 




rdlllliy ^JUbll IUIUI 1 1 

cluster 1 


Fr\ r\ a in rn 1 
rui vvui u 

Reverse 


1 6S 
16S 


AAAf^AAr.ATTA ATArrr^rATA 

AAAOOAAOA 1 1 AM 1 Av_^VJ^_A 1 A 

TCmXTAATCTCTACGCA 


1 59 
696 


538 


IVIUU 1 1 IcU IIUIII 

Rinttila et al. [37] 
Hung et al. [38] 


6 


Species Clostridium 
perfringens 


Forward 
Reverse 


16S 
16S 


TGAAAGATGGCATCATCATTCAAC 
GGTACCGTCATTATOTCCCCAAA 


183 
440 


258 


Skanseng et al. [39] 
Skanseng et al. [39] 


6 


Family Clostridium 
cluster IV 


Forward 


16S 


ACAATAAGTAATCCACCTGG 


866 


298 


Modified from 
Ramirez-Farias et al. [40] 


9 




Reverse 


16S 


CTTCCTCCCl 1 1 ICnICAA 

v_ 1 1 v_ v_ 1 v_ v_ vJ 1 1 1 1 \J 1 


1 163 




MatM iki pf nl T411 

iviaLjUr\i ct ui. L ' 1 J 




Famik/ r~lpictriHii im 
rdlllliy ^JUbll IUIUI 1 1 


Fr^r\Ainrn l 
rui VVUI U 


1 6S 


CCCT A CCTC A CTA ACAACC 

Luu 1 A^.^. 1 UnL 1 AAOAAO^_ 


482 


41 3 


Rinttilpi pt nl T371 

nl 1 1 LLI Id Ct Ul. ID / J 


13 


cluster XIV 


Qompvzp 
ncVcl ic 


1 ^ 


CTTTC A CTTTC ATTCTTCCC A A 

^.111 unu 1 1 1 \-.r\ 1 1 ^ 1 1 \3^.\3r\r\ 






K/latci \V\ Pt nl TA91 
IVIdLbUKI ct Ul. L^^-J 




Phvli im Rrar-torr^iHotoc 
rliyiUIII DdLLcl UlUcLcb 


Fni r\ a in rn 1 
rui VVUI U 


1 ^ 


CCCCA\AfTYATTCCCTTTA A ACCC 

^uUnVV 1 i A 1 1 \J\J\J 1 1 1 nnnUUU 




41 4 


K/li"ihlinn Pt nl T331 
IVIUIIIIIiy ct Ul. \_DD\ 


1 




Reverse 


16S 


GGTAAGGTOCTCGCGTA 


967 




Muhling et al. [33] 




otrl lUb DULLtzIUIUcb 


Fni r\ a in rn 1 
rui VVUI U 


1 ^ 


A A CCTCCCCC A C ATTCC 


DUZ 


3DD 

DUU 


N/lanv Pt nl TA31 
i Via 1 IZ ct Ul. \f-rD\ 


1 




Reverse 


16S 


GAGCCGCAAACmCACAA 


601 




Franks et al. [44] 




Phylum Actinobacteria 


Forward 


16S 


GCGKCCTATCAGCTO^ 


235 


333 


Modified from Glockner 
et al. [45] first published. 

1 1 Icy Iclcl LU Olid 11 L'UJ 


2 




Qompvzp 
ncVcl oc 


1 ^ 


rrr^rrTArnAnrYrTTTArnr 


DO/ 




Thic. cti irl\/ riiKtPil\A/9 T3^1 

II lib bLUUy, ^.lUbLdlVVZ \_D\J\ 




Family Bifidobacteriaceae 


Forward 


16S 


/ — i~/^~v i / — • / — AAA n~ / — ' / — rr~ f 

CTCCTGGAAACGGGTGG 


152 


442 


Matsuki et al. [42] 


2 




Reverse 


16S 


CmCACACCRGACGCG 


593 




Delroisse et al. [47] 




Class (3- and 
y-proteobacteria 


Forward 
Reverse 


23S 
23S 


GTATAATGGGTCAGCGAC 
CAGCATCGCAOTCTGA 


569 
1241 


673 


This study, Linux ARB [13] 
This study, Linux ARB [13] 


8 


Family Enterobacteriacea 


Forward 


16S 


CGTCGCAAGMMCAAAGAG 


182 


333 


Modified from 
Fried rich et al. [48] 


8 




Reverse 


16S 


TOCCGCGGCTGCTGGCAC 


514 




Modified from 
Palmer et al. [49] 




Species Escherichia coli 


Forward 


16S 


GTOATACCmGCTCATOA 


461 


320 


Malinen et al. [50] 


8 




Reverse 


16S 


ACCAGGGTATCTAATCCTG^ 


780 




Malinen et al. [50] 
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Table 1 Target-specific primer sequences constituting the Gut Micro biotassay in order to target main bacterial phyla 
and bacterial groups and species of high interest in the mammalian intestine (Continued) 



Class s-proteobacteria 


Forward 


16S 


TGGTGTAGGGGTAAAATCCG 


680 


286 


Bui etal. [51] 


5 




Reverse 


16S 


AGGTAAGGTOTOGYGTATC 


965 




This study, Primrose [14] 




Class 6-proteobacteria 


Forward 


16S 


GGTGTAGGAGTGAARTCCGT 


681 


534 


This study, Primrose [14] 


3 




Reverse 


16S 


TACGTGTGTAGCCCTRGRC 


1214 




This study, Primrose [14] 




Phylum Fusobacteria 


Forward 


16S 


GATCCAGCAATOTGTGTGC 


387 


292 


Sanguin et al. [52] 


10 




Reverse 


16S 


CGAAmCACCTCTACAOTGT 


678 




Walter et al. [53] 




Phylum Verrucomicrobia 


Forward 


16S 


GAATCTCGGTGTAGCA 


673 


551 


Modified from Ranjan [54] 


15 




Reverse 


16S 


GGCATOTAGTACGTGTGCA 


1223 




This study, Primrose [14] 




Phylum Spirochaetes 


Forward 


16S 


GTmAAGCATGCAAGTC 


45 


294 


Choi et al. [55] 


4 




Reverse 


16S 


TGCTGCCTCCCGTAGGAG 


338 




This study, ClustalW2 [36] 




Domain Archaea 


Forward 


16S 


CAGCMGCCGCGGTAAWC 


518 


440 


Giovannoni et al. [56] 


12 




Reverse 


16S 


YCCGGCGTOAMTCCAA^ 


957 




Delong [57] 





Reference bacteria listed are target for the respective primer sets. Numbers represent: 1: Bacteroides fragilis (DJF_B083(EU 728706)), 2: Bifidobacterium pseudolongum globosum 
(DMS 20092), 3: Bilophila wadsworthia (ATCC 49260), 4: Brachyspira pilosicoli (isolated from the intestine of slaughter pig at DTU-VET, 28-02-2000), 5: Campylobacter fetus 
(ATTC 10852), 6: Clostridium perfringens (NCTC 10240), 7: Enterococcus faecalis (ATCC 2921 2), 8: Escherichia coli (971 1 108-2), 9: Faecalibacterium prausnitzii (DSM 1 7677), 
10: Fusobacterium Necrophorum (ATCC 25286), 1 1: Lactobacillus sakei (DSM 2001 7), 12: Methanocorpusculum labreanum (DSM 4855), 13: Roseburia sp. (DJF_VR77 
(EU728794)), 14: Streptococcus suis (NCTC 10446), and 75: Verrucomicrobium spinosum (DSM 41 36). Nucleotide explanation: Y, CfT ; R, A/G ; W, A/T ; K, GfT ; M, A/C. 
*To improve coverage of the Genus Enterococcus, two reverse primers were mixed in equal concentrations. 



(NanoDrop Technologies Inc., Wilmington, Germany) 
spectrophotometer (Additional file 2). DNA was stored 
at -20°C until needed. 

Verifying the Gut Microbiotassay on the Access Array 
48.48 Tenfold serial dilutions ranging from 0.5 pg/|il 
to 50 ng/(il were made from DNA extracted from 15 
reference bacteria. The AA48.48 was processed follow- 
ing the Access Array System™ User Guide' [8]. In short: 
primer sets were mixed in equal concentrations and di- 
luted to 4 \iM with 20 x Access Array Loading Reagent 
(Fluidigm, South San Francisco, CA, USA) and nuclease- 
free water (Ambion Inc., Austin, USA). Master Mix was 
prepared as described in the instructions, and sample 
solutions were produced from 4 \A Master Mix and 1 \A 
DNA. The AA48.48 was primed in a pre-PCR' IFC control- 
ler AX (Fluidigm), before it was processed in the Biomark 
(Fluidigm) using the Fluidigm AA 48 x 48 Standard 
vl'-protocol listed in Table 2. Subsequently the ampli- 
cons were harvested in a 'post-PCR' IFC controller AX 
(Fluidigm). 

Harvested amplicons were measured on the Agilent DNA 
1000 chip (Agilent Technologies, Waldbronn, Germany) 
with the Agilent 2100 Bioanalyzer (Agilent Technologies) 
to verify the specificity of the Gut Microbiotassay. This was 
assessed from the size and number of amplicons generated 
by the primer sets listed in Table 1, with the reference 
bacteria as targets. 

Three primer concentrations (1, 2 and 4 |iM), and three 
primer:454BL-ratios (1:0.2, 1:0.4 and 1:0.8 \iM) were tested, 
before settling on the final protocol. 



Table 2 The Fluidigm 'AA 48 x 48 Standard vl' PCR 
thermal protocol [8] 



PCR stages 


Temperature 


Number 
of Cycles 


Thermal mixing and 
hot start phase 


50°C 2 minutes 
70°C 20 minutes 
95°C 10 minutes 
95°C 15 seconds 


1 


PCR cycle 


60°C 30 seconds 
72°C 1 minute 


10 


COt cycle 


95°C 15 seconds 
80°C 30 seconds 
60°C 30 seconds 


2 


PCR cycle 


72°C 1 minute 
95°C 15 seconds 
60°C 30 seconds 
72°C 1 minute 


8 


COt cycle 


95°C 15 seconds 
80°C 30 seconds 
60°C 30 seconds 
72°C 1 minute 


2 


PCR cycle 


95°C 15 seconds 
60°C 30 seconds 
72°C 1 minute 


8 


COt cycle 


95°C 15 seconds 
80°C 30 seconds 
60°C 30 seconds 


5 


Extension 


72°C 1 minute 


1 
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Testing the Gut Microbiotassay on complex samples 
Samples and sampling 

To explore the efficiency range of the Gut Microbiotassay, 
and test its sensitivities under varying circumstances, the 
study included 12 three-day-old piglets from conventional 
pig-farms (Danish landrace); seven with and five without 
clinical diarrhoea (Additional file 3: Table S2). None of the 
animals had received any treatment. Piglets were sacrificed 
and luminal content from the small and large intestine 
were obtained immediately after and frozen at -80°C All 
handling of the animals was in accordance with regulations 
from The Danish Centre for Animal Welfare. 

DNA extraction 

100 mg intestinal content was suspended in 900 ul PBS 
and bead-beated in 2 ml microcentrifuge tubes con- 
taining a 5 mm steel bead (Qiagen, Hilden, Germany) 
at 15.0 hertz for 2.3 min (Tissuelyser II, Qiagen). Tubes 
were centrifuged at 10,000 x g 90 s, 20°C, and 350 ul of 
supernatant were transferred to new tubes. These were 
placed in the QIAsymphony SP (Qiagen) for DNA ex- 
traction using the QIAsymphony Virus/Bacteria Mini 
Kit (Qiagen) with the protocol 'Pathogen complex 200' 
(Qiagen), elution volume: 60 ul. DNA was measured as 
previously mentioned and stored at -20°C until further 
processing (Additional file 3: Table S2). 

Analysing complex samples on the AA48.48 with the 
Gut Microbiotassay 

To test the amplification efficiency on DNA extracted from 
intestinal content sample-calibration-curves' were con- 
structed from tenfold serial dilutions (0.25 pg/ul - 25 ng/ul). 
These were made from pooled sample material of 1 ul 
extracted DNA from each of the samples listed in 
Additional file 3: Table S2. 'Control-calibration-curves were 
constructed from the IPC, and reference-calibration-curves' 
were generated by pooling all reference bacteria in equal 
concentrations. 

DNA extracted from the pig samples included in the 
study was diluted to 50 ng/ul with nuclease-free water 
(Ambion). To test how the addition of a 454BL affected 
the Cq values, an AA48.48 with a 454BL and one without a 
454BL were run under the exact same conditions. With 
a 454BL, each sample mixture consisted of 3 ul Master 
Mix, 1 ul sample (DNA, 50 ng/ul), and 1 ul Access Array 
Barcode Library for the 454 GS FLX Titanium Sequencer 
(Fluidigm), 2 uM. Afterwards the amplicons were harvested 
and stored at -20°C until needed. 

Amplicon preparation for NGS 

To normalize the harvested sample amplicons for NGS, 
they were measured as described earlier and pooled in 
equal concentrations. This resulted in a total volume of 
144.27 ul of 16 ng/ul which was purified to remove any 



PCR by-products. First the volume was reduced to 15 ul 
by extracting the DNA with phenol chloroform using a 
standard procedure [15]. Next, the extracted DNA was 
run in a 0.7% Seakem® LE Agorose gel (Lonza Rockland, 
Rockland, ME, USA) for 86 min, 90 V, and incubated for 
30 min in 0.0004% ethidium bromide. Bands were visual- 
ized with UV-radiation using the Bio-Rad Universal 
hood II (Segrate, Milan, Italy) and bands in the size range 
200-900 bp were excised, equalling expected amplicon 
sizes. DNA was extracted from the excised gel using the 
Qiaquick Gel Extraction Kit (Qiagen) in accordance with 
the kit manual. 

The final pool of purified DNA from the 12 piglets to 
be sequenced was 723.25 ng (260/280-nm ratio: 1.96). 
This was run on a quarter PicoTiterPlate™ on a 454 GS 
FLX Titanium Sequencer (Roche) by LGC Genomics 
(GmbH, Berlin, Germany). 

Data analysis 

From raw Cq values to relative quantitative data Cq 

values generated from the AA48.48 without a 454BL 
added were exported from the 'Fluidigm Real-Time PCR 
Analysis software', version 3.0.2 (Fluidigm), to Excel. Cq 
values were corrected to the IPC included in all runs, 
and those exceeding primer specific cut-off values, deter- 
mined from the verification step of the Gut Microbiotassay 
(Additional file 4: Table S3), were excluded. Relative quan- 
tification was calculated from the mean of the technical 
replicates using the Livak- method [16]. This method was 
chosen based on the theory that the total amount of 
bacteria, targeted by the primer set domain Bacteria B, 
constituted 100% of the microbiota in the individual 
gut section at all times. Hence, the Cq values of all pri- 
mer sets for each sample were normalised against the 
Cq value of their respective domain Bacteria B primer set: 
R(primersetx) = 2 Cq ( d — B)~ C( iij>rimer set x) # To compare total 
bacteria detected with the domain Bacteria B primer set, 
these Cq values were related to the final number of 
thermal cycles run, 35: R (domain Bacteria b) = 2 35 - Cq (^Bac t e r iaB) t 
Normalization was further done to total mean of all primer 
sets for each sample with similar results (data not shown). 

454- sequencing data Sequence data, available at NCBI 
Sequence Read Archive under Accession SRA061551, 
was analysed using BION, a yet unpublished open source 
program. For more information on the BION software, its 
functions, and the main statistics for the raw results see 
Additional file 5. In short: the sequence dataset was con- 
verted to FASTQ, split according to sample barcodes and 
primer sequences, and trimmed in the same process. Next, 
sequences were cleaned using a cut-off for minimum qual- 
ity of 96%, and a minimum sequence length of 200 bp. The 
remaining sequences were clustered using UCLUST, based 
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on a minimum seed similarity of 99.5%. Query sequences 
were compared to Greengenes Gene Database [17] using 
the k-mer matching program Simrank2, an improved 
version of Simrank [18]. Simrank2 returns the n % best 
similarities, no matter how low they are. Also, it only 
produces 8-mers from regions above a given quality, and 
skips sequences with too few 8-mers. This last feature 
helps improve the data quality. Simrank2 was set to return 
the best 1% similarities with a similarity cut-off of 50%. 
Taxonomy was generated by transforming percentages for 
the Greengenes OTUs to scores, read densities, the sum 
of which was 1. This was done in a weighted manner, 
so that OTUs of the highest similarity scored a high num- 
ber and vice versa, the sum of which was the original 
reads. Only phylotypes with a primer-specific read density 
of > 1% were included in the statistical analysis. 

Statistics Initially, all data were transformed with the 
natural logarithm. Primer sets with more than half of the 
data missing were removed from the dataset. If primer 
sets in taxonomic lineage showed pair-wise correlations 
above 0.99, only the primer set of the highest taxonomic 
level in the lineage was retained for analysis, and any con- 
clusion drawn from this also accounted for the excluded 
sub-level primer sets. To include the information from 
remaining data-deficient primer sets, for which the fraction 
of missing data was low and never above 0.35, the EM- 
algorithm [19] was used to substitute missing values with 
imputed ones, by applying a multivariate Gaussian model. 
Each primer set was allowed to depend freely on the others 
and also for dependency on gut section, diarrhoeic status, 
and interaction between these. Model fit of the multivariate 
model for the primer sets after imputation was assessed by 
transforming the model residuals with the inverse of the 
square root of the estimated covariance matrix between the 
primer sets, and applying standard model control to these 
standardized residuals. These analyses were consistent with 
standard model behaviour. 

Effect of gut section and diarrhoeic status was tested with 
the Likelihood Ratio Method using the Wilks test [20]. 
Sequence data were compared for effect of diarrhoeic 
status using the non-parametric Wilcoxon test [21]. 

P-values < 0.05 were considered statistically significant. 

Results 

Designing the Gut Microbiotassay 

The Gut Microbiotassay was constructed from 49 primers 
constituting 24 primers sets (Table 1) in order to target the 
main bacterial Phyla (Actinobacteria, Bacteroides, Firmicutes, 
Fusobacteria, Proteobacteria and Verrucomicrobia) and 
some of the highly important bacterial groups and species 
reported in the mammalian intestinal tract [9]. 12 primers 
were designed de novo, and 37 primers were from pub- 
lished literature, of which ten were modified. Five of 



the 24 primer sets were unmodified pairs from published 
literature, whereas the rest were combinations composed in 
this study. As it was not possible to find or design specific 
primer sets targeting the 16S rRNA gene for class (3- and 
y-proteobacteria, and genus Lactobacillus, these were 
designed to target the 23S rRNA gene instead. 

Sensitivity and specificity of the Gut Microbiotassay 

Tenfold serial dilutions of DNA extracted from 15 ref- 
erence bacteria were used to evaluate the specificity 
and sensitivity of the Gut Microbiotassay (Figure 1). 

The specificity for each primer set was assessed from the 
Cq values obtained for their respective target bacteria in 
addition to any cross reaction. For 13 of the 24 primer sets, 
specific positive reactions were registered on various taxo- 
nomic levels in agreement with the reference bacteria listed 
for each primer set in Table 1. The remaining primer sets 
showed different degrees of cross reaction (Figure 1). 
On average, the lowest unspecific and the highest spe- 
cific Cq of same concentration differed by 9 Cq values, 
equalling less than 1% cross reaction. Apparently, the 
primer set phylum Firmicutes did not have a complete 
coverage, as it only amplified few bacteria from the 
Clostridia clusters. The highest specific Cq value was 
determined for each primer set from the representative 
target bacterial species and used as cut-off value in the 
data analysis (Additional file 4: Table S3). 

The dynamic range for the specific primer sets spanned 
from 50 ng to 50 pg DNA/|il for two primer sets, to 5 pg 
DNA/(il for 15 primer sets, and for seven primer sets down 
to 0.5 pg DNA/(il (lowest concentration tested), Table 3. 
Linear regression of log-concentration versus Cq values 
for the different primer sets demonstrated r 2 from 0.943 
to 0.999 (SD 0.0074). Gel pictures and electropherograms 
confirmed specific amplification by demonstrating the 
expected number of amplicons corresponding to the num- 
ber of taxonomic levels represented by the specific primer 
sets. Additionally, amplicon sizes (bp) were comparable to 
those listed in Table 1. 

Fluidigms recommendation on DNA sample concentration 
is 25-50 ng/(il [8]. When testing the amplification efficiency 
on 23 calibration curves (ranges: 0.25 pg/|il - 50 ng/|il), 
there was no PCR inhibition for either of the primer sets 
with DNA of 25 or 50 ng/ul 

Addition of a 454 Barcode Library 

Including a 454BL with the reactions did not affect the 
Cq values drastically. The biggest deviation was less than 1 
Cq compared to the respective Cq values without a 454BL 
added, and 90 percent of the corresponding values from 
the two data sets were correlated by r 2 > 0.90. Uncorrelated 
data were mainly seen if no or very few target bacteria were 
present in the reaction. In such cases the 454BL occasion- 
ally caused some interaction with itself or/and the present 
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Tenfold serial dilution of DNA extracted from reference bacteria (50 ng/u.1 - 0.50 pg/ul) 



Domain Bacteria A- 
Domain Bacteria B- 
Phylum Firmicutes- 
Class Bacilli- 
Genus Enterococcus - 
Genus Streptococcus - 
Genus Lactobacillus - 
Family Clostridium cluster I - 
Species Clostridium perfringens - 
Family Clostridium cluster IV - 
Family Clostridium cluster XIV - 
Phylum Bacteroidetes - 
Genus Bacteroides - 
Phylum Actinobacteria - 
Family Bifidobacteriaceae- 
Phylum Fusobacteria - 
Class B- and y-proteobacteria - 
Family Enterobacteriaceae- 
Species Escherichia coli - 
Class 6-proteobacteria- 
Class e-proteobacteria - 
Phylum Verrucomicrobia - 
Phylum Spirochaetes- 
Domain Archaea - 



■ ■ 
■ 

■ | ■ 




Cq value 



■ 


<1 no 




<10,62 




11,22 




12,44 




13,65 




14,86 




16,08 




17,29 




18,50 




19,72 




20,93 




22,15 


5 


23.36 


■ 


24,57 


■ 


25,79 


■ 


27,00 


■ 


28,21 


■ 


29,43 


■ 


30,64 


■ 


31,86 


■ 


> 32,46 


■ 


> 35,00 



Figure 1 Heatmap generated by the Fluidigm Real-Time PCR Analysis program from raw Cq data. At the top horizontal: Tenfold serial 
dilutions of DNA (ranging from 50 ng/ul down to 0.5 pg/ul) extracted from 15 different reference bacteria. Raised numbers in parentheses 
represent the respective numbers given to each reference bacterium in Table 1. These were used to test the specificity and sensitivity of the 
primer sets included in the Gut Microbiotassay, listed vertically on the left. Primers were run in duplicates. On the right: a colour scale depicting 
the respective Cq values (the software uses commas instead of points as decimal separator for Cq values). 



primer set, creating a signal that obscured the low or 
missing signal from the target bacteria. Because of this, 
samples were run twice: with and without a 454BL added. 
Cq values without a 454BL were used for data analysis, 
whereas amplicons with the 454BL incorporated were 
available for 454-sequencing. 

The Gut Microbiotassay on complex samples 

Piglets of different diarrhoeic status were sacrificed and 
luminal content was collected from the small and the large 
intestine. This resulted in a total of 23 samples divided into 
the following groups: small intestine without diarrhoea (S-) 
n = 5, small intestine with diarrhoea (S+) n = 7, large in- 
testine without diarrhoea (L-) n = 4, and large intestine 
with diarrhoea (L+) n = 7. 

AA48.48 general findings 

The following primer sets were missing more than 
half of the data, and were consequently removed from 
the statistical data analysis: family Clostridium cluster 
IV, class e-proteobacteria and 5-proteobacteria, phyla 
Verrucomicrobia, and domain Archeae. Primer sets 
showing pair-wise correlation above 0.99, were: class 



(3- and y- proteobacteria with family Enterobacteriacea, 
and species E. coli; and phylum Bacteroidetes with genus 
Bacteroides. Of the general bacteria primers only domain 
Bacteria B was tested as this was the primer set used for 
normalization. This left 15 primer sets for data analysis. 
The multivariate test for effect of gut section and diarrhoeic 
status, with effects allowed to differ with the values of gut 
section, revealed statistically significant effects (S/L p = 0.01, 
S+/S-p =0.002, and L+/L-p =0.006, respectively). 
However, effects were not limited to the primer sets that 
showed significant differences in the Gut Microbiotassay. If 
these were excluded, and the data analysis repeated, this 
still resulted in a significant effect from diarrhoeic status in 
the large intestine (L+/ L - p = 0.01), while it was borderline 
insignificant in the small intestine (S+/ S -p = 0.056). 

Table 4 lists the estimated mean percentages and 95% 
confidence interval for each primer set included in the 
data analysis relative to total bacteria. Diarrhoeic gut 
sections contained fewer bacteria from phylum Firmicutes 
(genus Streptococcus and unclassified (p < 0.05)), but 
a higher fraction of genus Enterococcus (small intestine 
p = 0.04). There was a highly significant depletion of mem- 
bers from genus Streptococcus in the diseased compared to 
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Table 3 The primer efficiency, revalue, dynamic range 
and limit of detection of the Gut Microbiotassay 



Primer set 


Efficiency, 

% 


r 2 


Dynamic 
range 


Limit of 
detection, 
ng/|jl 


Domain Bacteria B V4-V5 


75-106 


>0.968 


4-6-fold 


0.05-0.0005 


Phylum Firmicutes* 


86-97 


>0.994 


3-5-fold 


0.05-0.005 


Phylum Bacilli 


99-101 


>0.943 


4-5-fold 


0.05-0.005 


Genus Enterococcus 


87 


>0.998 


5-fold 


0.005 


Genus Lactobacillus 


112 


>0.992 


5-fold 


0.005 


Genus Streptococcus 


102 


>0.999 


5-fold 


0.005 


Family Clostridium cluster 1 


96 


>0.995 


4-fold 


0.05 


Species Clostridium perfringens 


87 


>0.999 


5-fold 


0.005 


Family Clostridium cluster IV 


103 


>0.997 


5-fold 


0.005 


Family Clostridium cluster XIV 


96 


>0.999 


4-fold 


0.05 


Phylum Bacteroidetes 


91 


>0.991 


6-fold 


0.0005 


Genus Bacteroides 


85 


>0.997 


5-fold 


0.005 


Phylum Actinobacteria 


92 


>0.988 


6-fold 


0.0005 


Family Bifidobacteriaceae 


77 


>0.995 


5-fold 


0.005 


Class (3- and y-proteobacteria 


87 


>0.999 


5-fold 


0.005 


Family Enterobacteriacea 


84 


>0.996 


5-fold 


0.005 


Species Escherichia coli 


91 


>0.999 


6-fold 


0.0005 


Class e-proteobacteria 


86 


>0.995 


6-fold 


0.0005 


Class 6-proteobacteria 


91 


>0.995 


5-fold 


0.005 


Phylum Fusobacteria 


85 


>0.997 


4-fold 


0.005 


Phylum Verrucomicrobia 


107 


>0.990 


6-fold 


0.0005 


Phylum Spirochaetes 


80 


>0.996 


5-fold 


0.005 


Domain Archaea 


95 


>0.996 


5-fold 


0.005 



*Data from the reference bacteria Roseburia sp. and F. Prausnitzii has been 
excluded in primer set Phylum Firmicutes. 

The Gut Microbiotassay was tested against tenfold dilution series of DNA extracted 
from individual reference bacteria without a 454 Barcode Library added. 



the healthy intestines. The same tendency was seen for 
Clostridium cluster I (small intestine p = 0.02). The diar- 
rhoeic small intestine harboured fewer members from class 
|3- and y-proteobacteria than the healthy one, while this 
was reversed for the large intestine (p < 0.05). Generally 
healthy piglets had a gut microbiota dominated by 
Gram-positive bacteria, which was partly displaced by 
Gram-negative bacteria in diarrhoeic piglets. 

454-Sequencing results 

Sequencing barcoded sample amplicons resulted in 275,133 
unprocessed consensuses, which dropped to 164,055 after 
quality trimming. Amplicons generated by the primer sets 
domain Bacteria A and B encompassed 16S rRNA gene 
sequences from both Gram-positive and Gram-negative 
bacteria genera plus from unclassified bacteria. Good 
congruency was found between the primer sets of the 



Gut Microbiotassay and the 454-sequencing results gener- 
ated from their respective amplicons (Figure 2). 

Comparing significant findings from AA48.48 with 
454-sequencing data 

Based on the Gut Microbiotassay data for the individual 
bacterial groups, four groups showed significant differ- 
ences: domain Bacteria B, genus Streptococcus, phylum 
Actinobacteria, and class |3- and y-proteobacteria, Table 4. 
In the small intestine, diarrhoea was associated with a 
significantly reduced bacterial load {p = 0.003). Se- 
quence data generated by the domain Bacteria B pri- 
mer set indicated that diarrhoea was associated with a 
decreased number of members from genus Streptococcus, 
including S. alactolyticus in both the small intestine 
(p = 0.0061, and p = 0.0099, respectively) and the large 
intestine (p = 0.0061 and p = 0.013). When comparing 
the results for genus Streptococcus, there were sig- 
nificantly fewer Streptococci in the diseased intestines 
(S+/ S-p = 0.0003, L+/ L-p = 0.00002). The sequences 
generated by genus Streptococcus confirmed the previous 
results (S+/ S -p = 0.0061, L+/ L-p = 0.012). These were 
classified to: genus Streptococcus, species S. hyointestinalis, 
S. suis, and S. alactolyticus. However, at species level, only 
S. alactolyticus, showed coherently significant results with 
genus Streptococcus (S+/ S -p = 0.047, L+/ L -p = 0.012). 
Phylum Actinobacteria was more abundant in the diar- 
rhoeic compared to the healthy small intestine (p = 0.034), 
however, this was not reflected in the sequencing data 
for the small intestine. The diarrhoeic large intestine 
possessed significantly more members from the class (3- 
and y-proteobacteria, including family Enterobacteriaceae 
and species E. coli, than the healthy one (p = 0.029). 
Nonetheless, for this gut section no significant differ- 
ences showed up in the sequencing data from the family 
Enterobacteriaceae primer set. 

Specificity of the Gut Microbiotassay 

After sequencing the amplicons generated by the Gut 
Microbiotassay the specificity of the primer sets was 
revaluated (Figure 2), except for the primer sets targeting 
the 23S rRNA gene (genus Lactobacillus and class (3- and y- 
proteobacteria) as the freeware used to analyze the sequence 
data is currently based on 16S rRNA gene databases. 

16 of the 24 primer sets only generated sequences 
from their intended target group, four showed cross re- 
action with a single group, and two primer sets cross 
reacted with more than two other taxonomic groups. 
The primer set phylum Firmicutes did not reveal any 
Clostridia, as predicted earlier. In conclusion, the sequen- 
cing data confirmed the specificity of the primer sets found 
in the verification process of the Gut Microbiotassay 
when tested against DNA extracted from pure-cultured 
reference bacteria. 
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Table 4 Estimated mean primer set values from the AA48.48 and the corresponding significant 454 GS FLX Titanium 
sequencing results 



Primer set 


Estimated mean percentages relative to total bacteria 


P-values 


S- 


S+ 


L- 


L+ 


S+/ S- 


L+/ L- 


Domain Bacteria B 1 


24.96 


20.32 


25.13 


25.46 


0.003 


0.857 




[22.39;27.73] 


[16.83:24.32] 


[22.26:28.27] 


[23.23:27.84] 






Genus Streptococcus 










0.0061 


0.0061 


Species S. alactolyticus 










0.0099 


0.0134 


Phylum Firmicutes 


36.08 


27.59 


30.59 


17.79 


0.649 


0.278 




[17.64;65.91] 


[7.50;72.58] 


[13.64:59.56] 


[9.79;29.83] 






Genus Streptococcus 










0.0242 


0.0102 


Unclassified 










0.0462 


0.0176 


Class Bacilli 


14.01 


15.92 


8.30 


12.04 


0.473 


0.52 




[4.85;32.09] 


[2.22;57.40] 


[2.50;20.64] 


[4.99;24.63] 






Genus Streptococcus 










0.0287 


0.0061 


Unclassified 










0.1516 


0.029 


Genus Enterococcus 


0.036 


0.39 


0.04 


0.11 


0.208 


0.34 




[0.0038;0.14] 


[0.0049:2.51] 


[0.0032:0.18] 


[0.018:0.39] 






Genus Enterococcus 










0.0424 


0.9273 


Genus Streptococcus 


4.04 


1.00 


5.10 


0.56 


0.0003 


0.00002 




[1.84;7.73] 


[0.24;2.83] 


[2.10;10.46] 


[0.29:0.98] 






Genus Streptococcus 










0.0061 


0.0121 


Species S. alactolyticus 










0.0467 


0.0121 


Genus Lactobacillus 


33.97 


53.70 


33.27 


27.55 


0.643 


0.846 




[12.03;76.74] 


[7.82:190.09] 


[10.27:81.52] 


[11.62:55.67] 






Family Clostridium cluster 1 


0.65 


0.17 


0.32 


0.079 


0.231 


0.139 




[0.16;1.81] 


[0.012;0.76] 


[0.064:0.97] 


[0.025:0.19] 






Genus Clostridium 










0.0171 


0.2558 


Species Clostridium perfringens 


3.03 


2.82 


2.33 


0.63 


0.81 


0.219 




[0.62;9.24] 


[0.14;14.31] 


[0.38;7.86] 


[0.17;1.66] 






Family Clostridium cluster XIV 


0.89 


16.84 


10.82 


5.02 


0.818 


0.672 




[0.098:3.53] 


[0.22:106.79] 


[0.86;47.63] 


[0.82;16.95] 






Phylum Bacteroidetes and Genus Bacteroides 


2.25 


76.17 


22.05 


5.10 


0.924 


0.567 




[0.082:12.19] 


[0.091:505.28] 


[0.48:130.66] 


[0.34:23.65] 






Species Bacteroides uniformis (generated from the primer set Genus Bacteroides) 






0.0447 


0.0985 


Phylum Actinobacteria 


0.32 


3.92 


1.02 


0.32 


0.034 


0.068 




[0.13:0.68] 


[0.70;12.70] 


[0.35;2.32] 


[0.15:0.61] 






Genus Actinomyces 










1 


0.0363 


Species Actinomyces hyovaginalis 










0.1042 


0.0179 


Species Bifidobacterium pseudolongum 










1 


0.0368 


Unclassified 










0.072 


0.0104 


Family Bifidobacteriaceae 


0.0064 


0.011 


0.0063 


0.012 


0.399 


0.24 




[0.0026:0.013] 


[0.0022:0.036] 


[0.0023:0.014] 


[0.0056:0.022] 






Species Bifidobacterium pseudolongum 










1 


0.0424 


Phylum Fusobacteria 


1.59 


32.60 


6.28 


9.71 


0.413 


0.604 




[0.18:6.17] 


[0.48:203.87] 


[0.54;27.06] 


[1.66:32.13] 
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Table 4 Estimated mean primer set values from the AA48.48 and the corresponding significant 454 GS FLX Titanium 
sequencing results (Continued) 



Class (3- and v-proteobacteria; Family 2 96 0.54 
Enterobacteriaceae and Species Escherichia coli [0 83;7 69] [0.049,2.28] 


0.61 
[0.14:1.72] 


3.17 
[1.10;7.25] 


0.6 


0.029 


Family Enterobacteriaceae (generated from the primer set Family Enterobacteriaceae) 






0.0424 


1 


Species Escherichia fergusonii (generated from the primer set Family Enterobacteriaceae) 






0.0368 


0.7758 


Genus Shigella (generated from the primer set Family Enterobacteriaceae) 






0.0462 


0.9244 


Species Shigella flexneri (generated from the primer set Family Enterobacteriaceae) 






0.0413 


0.5626 


Phylum Spirochaetes 0 0059 0 010 


0.0091 


0.0056 


0.73 


0.081 


[0.004;0.0085] [0.005;0.019] 


[0.0058:0.014] 


[0.004:0.0076] 






Genus Streptococcus 






0.0169 


0.0127 


Unclassified 






0.0191 


0.0694 



domain Bacteria B is expressed in percentage relative to the final number of thermal cycles run, 35. 

Primer set values from the AA48.48 are expressed as percentages relative to total bacteria, and beneath (indented) are listed the significant 454-sequencing 
results for the respective primer sets. Numbers in bold indicate significant p-values. L: Large intestine, S: small intestine, +: with diarrhoea, -: without diarrhoea. 
[95% lower level; 95% upper level]. 



Discussion 

Most qPCR studies describing the gut microbiota typic- 
ally do this by using a general bacteria primer and a few 
group-specific ones [4,22]. In contrast, this study devel- 
oped the Gut Microbiotassay, an assay composed of 24 
primer systems, capable of screening the microbiota for 
the most common bacteria in the mammalian intestine 
[9,10] at various taxonomic levels. The Gut Microbiotassay 
was tested against representative reference bacteria, and 
next on complex intestinal samples from piglets of different 
diarrhoeic status. The sample amplicons were harvested 
and sequenced, and functioned as a proof of concept, 
evaluated the specificity of the Gut Microbiotassay by 
further elucidating the components of the gut microbiota. 

This approach offers an alternative to current molecular 
methods employed to characterize the gut microbiota such 
as phylogenetic microarrays [23], and NGS [24]. In contrast 
to phylogenetic microarrays, the AA48.48 is highly flexible 
because primer sets can readily be replaced to meet the 
needs of a current research study. In addition, the AA48.48 
outmatches the phylogenetic microarray on sample cap- 
acity, as well as sensitivity [25]. Also, no pre-amplification 
is needed when running the Gut Microbiotassay with the 
AA48.48, which reduces the workload, and also the risk of 
introducing technical variation. The effect of such technical 
variation can be reduced by normalization. In the present 
study each sample were normalized against the Cq value of 
their respective domain Bacteria B primer set. Impact of 
normalization against their respective domain Bacteria 
B was tested by performing a second normalization pro- 
cedure; at this point data was also normalized to total 
mean of all primer sets for each sample individually, 
with similar end results (data not shown). As the choice 
of normalization in the present study (Domain Bacteria B) 
or total mean did not have a great impact on final results 
we assume that the technical variation was low or that we 



managed to normalize efficiently using both methods. Fur- 
ther, data from domain Bacteria A and domain Bacteria B 
were highly correlated also pointing to domain Bacteria B 
as a reasonable reference primer efficiency of qPCR is of 
great importance and should ideally be within the range of 
85% to 110%. The majority of the primers used in the Gut 
Microbiotassay (21 out of the 23 primers included in 
Table 3) had and acceptable efficiency between 80-110%. 
However, the primer Doman bacteria A was found to dif- 
fer to much with regards to efficiency and dynamic range, 
thus data from this primer pair was only used to support 
data generated from the primer Doman bacteria B. Likewise 
primer set Phylum Firmicutes produced to high efficiency 
when tested on two reference bacteria Roseburia sp. and 
F. Prausnitzii. However Roseburia sp. and F. Prausnitzii are 
covered by several other well performing primers including 
Clostridium cluster XIV and Clostridium cluster IV respect- 
ively. Therefore the Gut Microbiotassay cannot be used as 
an absolute quantitative assay across the different primer 
sets. In order to make it a truly quantitative assay it will be 
necessary either to have a defined start sample regarding 
species composition or only to use specific primers on the 
array where high efficiency has been proven. 

The heatmap generated by the software 'Fluidigm 
Real-Time PCR Analysis' (Fluidigm Corporation) depicts 
the raw Cq values for each reaction. This makes it feas- 
ible to quickly evaluate and visually compare the bacter- 
ial profiles across a large number of samples. As sample 
amplicons are harvested individually following qPCR, 
it is possible to pinpoint which samples to sequence 
for further taxonomic information. Selective sequen- 
cing reduces costs compared to non-selective sequencing. 
Also, the dataset generated from sequencing the Gut 
Microbiotassay PCR amplicons produces a much more 
manageable dataset compared to metagenomic approaches. 
The Gut Microbiotassay provides a quantitative picture 
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Figure 2 Specificity of the Gut Microbiotassay. The figure shows the sequencing results and corresponding read density scores (> 1%) for the individual primer sets of the Gut Microbiotassay 
targeting the 16S rRNA gene. Green colours indicate sequences corresponding to the primer sets' target, whereas red colours represents unspecific amplification. The darker the colour, the higher 
the read density score. 
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of the distribution of the known gut microbiota repre- 
sented by the primers. Moreover, if combined with 
454-sequencing, it enables detection of bacteria with un- 
identified sequences [25]. A limitation of the assay when 
running on complex bacterial samples is that the primers 
will have different efficiencies and dynamic ranges due to 
imperfect matches with some of the target sequences. The 
Gut Microbiotassay has therefore the most value for analys- 
ing high-throughput quantification of the bacterial compos- 
ition in many samples or samples with defined biomarkers. 

The validation of the Gut Microbiotassay by sequencing 
the amplicons from intestinal content of two different gut 
sections from piglets of different diarrhoeic status demon- 
strated the potential to further elucidate the components 
of the gut microbiota. This study used the results from the 
Gut Microbiotassay to quantify the taxonomical groups, 
and NGS to access the bacterial constituents. 

Common intestinal bacteria in the neonatal piglet include 
members of: Clostridia, Streptococcaceae, Lactobacillaceae, 
Enterobacteriaceae, Fusobacteria and sometimes Bacter- 
oidetes [26,27]. These bacterial groups were also found 
in the gut microbiota of three-day-old piglets using the 
Gut Microbiotassay. The Gut Microbiotassay indicated 
that a healthy gut microbiota was dominated by Gram- 
positive bacteria, which were partly replaced by Gram- 
negative bacteria in the large intestine of diarrhoeic pig- 
lets. Robinson et al. [28] came to a similar conclusion in a 
study investigating the intestinal microbiota of pig colons 
experimentally induced with swine dysentery. Consistent 
significant findings from the Gut Microbiotassay and the 
454-sequencing results implied that diarrhoea was as- 
sociated with a depletion of members from the genus 
Streptococcus, and previous research has shown that 
Streptococci is an important member of a healthy gut 
microbiota [26,28,29]. A detailed review of the aetiology 
behind the piglet diarrhoea is beyond the scope of 
this paper and prevented by the limited number of 
piglets analysed, as the primary focus of this paper has 
been on the verification and application of the Gut 
Microbiotassay. 

Conclusions 

The Gut Microbiotassay offers affordable quantitative 
screening of the microbiota with the AA48.48. It has been 
thoroughly tested and strict criteria for data analysis have 
been outlined. It provides a high sample capacity, a wide 
dynamic range, and it facilitates selective 454-sequencing 
afterwards. Hence, it is timesaving and economical due to 
the easy library preparation, the low consumption of master 
mix, and the optional selective sequencing. These features 
make the Gut Microbiotassay a worthy high-throughput 
competitor to the current alternative methods used for 
investigating diverse ecosystems. 
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